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We calculate the spin-independent scattering cross section for direct detection that results from the 
electromagnetic polarizability of a composite scalar baryon dark matter candidate - “Stealth Dark Matter”, 
based on a dark SU(4) confining gauge theory. In the nonrelativistic limit, electromagnetic polarizability 
proceeds through a dimension-7 interaction leading to a very small scattering cross section for dark matter 
with weak-scale masses. This represents a lower bound on the scattering cross section for composite 
dark matter theories with electromagnetically charged constituents. We carry out lattice calculations of 
the polarizability for the lightest baryons in SU(3) and SU(4) gauge theories using the background field 
method on quenched configurations. We find the polarizabilities of SU(3) and SU(4) to be comparable 
(within about 50%) normalized to the baryon mass, which is suggestive for extensions to larger SU(N) 
groups. The resulting scattering cross sections with a xenon target are shown to be potentially detectable 
in the dark matter mass range of about 200-700 GeV, where the lower bound is from the existing LUX 
constraint while the upper bound is the coherent neutrino background. Significant uncertainties in the 
cross section remain due to the more complicated interaction of the polarizablity operator with nuclear 
stmcture, however the steep dependence on the dark matter mass, Ijrrfg, suggests the observable dark 
matter mass range is not appreciably modified. We briefly highlight collider searches for the mesons in 
the theory as well as the indirect astrophysical effects that may also provide excellent probes of stealth 
dark matter. 

PACS numbers: ll.I5.Ha, 12.60.-i, 95.35.+d 


Introduction - Despite remarkable advances in direct 
detection experiments |[T}01, a conclusive signal of nuclear 
interactions with dark matter (DM) remains elusive. These 
experiments, wbicb are sensitive to nucleus-DM cross sec¬ 
tions at or below 10“"'® cm^ per nucleon, bave already ex¬ 
cluded large classes of interactions and models, and are 
now actively probing Higgs boson exchange Em. 

Composite DM, which arises as a neutral bound state 
of a strongly-coupled gauge force (for early work, see ||3- 
mi has sparked multiple recent lattice calculations ifTOl - 
[ni. If its constituents are electromagnetically charged, 
the DM will interact with standard model (SM) nuclei via 
photon exchange, with the cross section suppressed by a 
momentum-dependent electromagnetic form factor. Ex¬ 
panding at small momentum transfer, one can obtain a se¬ 
ries of effective operators describing the interaction: the 
dimension-5 magnetic moment, dimension-6 charge ra¬ 
dius, and dimension-7 polarizability are the leading opera¬ 
tors ifTSti^ . Symmetry considerations can give models in 
which the first two operators are identically zero ll^l22l . 


Scattering due to electromagnetic polarizability remains, 
giving a lower bound on the direct detection cross section 
for a composite DM particle with charged constituents. 

We will focus here on a particular composite DM model, 
“stealth dark matter” ll23l . in which the DM is a scalar 
baryon composed of dark fermions that transform under an 
SU(iV£)) theory with even. In this model, electroweak 
symmetry breaking proceeds through the SM Higgs mech¬ 
anism. The dominant modes of interaction with the SM 
are the polarizability operator and direct Higgs boson ex¬ 
change. The latter was studied in some detail in ifT^ 1^ . 
placing bounds on the allowed dark matter coupling to the 
Higgs boson. In this work we study the polarizability, 
which unlike the Higgs interaction has no adjustable pa¬ 
rameters, but rather is completely determined by the strong 
dynamics once the gauge group and matter content are 
specified. 

SU(4) baryons and polarizabilities - A full consfruc- 
fion of fhe sfealfh DM model is given in ll2^ : here we 
briefly summarize fhe relevanf defails. The dark sector con- 
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sists of an unbroken SU(4) gauge theory, which contains 
bosonic baryonic bound states made up of four constituent 
fermions. The DM candidate itself is a scalar made up of 
two pairs of fermions which are degenerate in mass and 
carry equal but opposite electric charges of ±1/2. Hence, 
there is no magnetic moment or charge radius, leaving just 
the electromagnetic polarizability as the dominant interac¬ 
tion with photons. 

Previous estimates of the polarizability of a composite 
scalar have led to direct-detection cross sections on the or¬ 
der of 10“^® cm^ II 20 II , approaching the interaction strength 
at which background neutrinos are expected to contaminate 
the DM recoil signal. However, the estimates were based 
on semi-classical calculations of a strongly-coupled inter¬ 
action, and so have uncontrolled uncertainties. 

Additionally, due to how internal electric charges are 
correlated, the polarizability of bosonic 4-fermion baryons 
may differ appreciably from QCD-based estimates. In one 
limit where the internal constituents are uncorrelated, the 
polarizabilities are expected to be comparable. However, 
if alternate flavors tend to form pairs based on their Pauli 
statistics, the 4-fermion baryon polarizability would be 
derivative-suppressed compared to the 3-fermion baryon 
(i.e. two dipoles vs. one dipole and one charge). In order 
to quantify this effect, we perform lattice calculations for 
both the SU(3) and SU(4) baryon polarizabilities. 

Polarizability and Direct Detection — The electric po¬ 
larizability of the scalar baryonic composite DM field B 
with mass ms can be written as an effective operator of 
the form 

OF = CFB*BF^^^Fy^v, ( 1 ) 

where F^^ is the electromagnetic field strength tensor, 
Vf_i = (1,0,0,0) in the static limit, and Cf is the polariz¬ 
ability with mass dimension —3 in the nonrelativistic limit. 
Only the electric polarizability is considered since the mag¬ 
netic polarizability is expected to be suppressed Il24^ . This 
is a two-photon vertex, so that the scattering off of nuclei 
will involve a virtual photon loop. Because this operator 
is induced at a high scale (roughly the dark confinement 
scale Af) ~ ms), it is expected to generate other interac¬ 
tions with SM particles when the appropriate effective field 
theory matching and running down to the nuclear scale are 
carried out ll25] - l28l : in fact, an explicit treatment for the 
polarizability operator is given in Il29ll . Although the ef¬ 
fects of the additional induced operators are not negligible 
in general, we And that they are small compared to the un¬ 
certainties (particularly from nuclear physics) and so we 
will omit them. 

From the interaction shown above, the coherent DM- 
nucleus scattering cross section (per nucleon) is given by 

a^ucleoniZ, A) = ^ (I Cf/)? f ) , (2) 

where m„ is the nucleon mass, fins = mniTiB /± 
ms) is the reduced mass, (Z, A) are the atomic and mass 


numbers of the target nucleus, and the angular brackets rep¬ 
resent the momentum-averaged form factors for heavy DM 
candidates in a given experiment ll29]l . 

The primary source of systematic uncertainty is on 
the nuclear physics side of the calculation - evaluat¬ 
ing the non-perturbative nuclear matrix element, fp = 
{AlF^"Fpi,\A). Various attempts to perturbatively esti¬ 
mate this matrix element have been performed with varying 
levels of complexity 1129^311 . But, the matrix element also 
has nontrivial excited-state structures that likely require a 
fully non-perturbative treatment. This matrix element is 
similar to those needed for double-beta decay experiments, 
estimates for which have substantial variation ||32j[33l. Un¬ 
til a more accurate extraction of this matrix element is per¬ 
formed, we will use dimensional analysis arising from non¬ 
relativistic loop momenta counting, 

/p ~ 3z (3) 

where R = l.2A^/^ fm, as used in the double beta de¬ 
cay context, a is the Ane-structure constant, and Mp is a 
dimensionless parameter. With the factor of 3 in Eq. 0, 
our expression approximately matches ll29l |30l for heavy 
nuclei when Mp ~ 1. To allow for an order of mag¬ 
nitude uncertainty in the nuclear matrix element, we take 
1/3 < Mp < 3, although a detailed nuclear structure 
extraction would be needed for a more precise estimate. 

Background field method - In order to extract the elec¬ 
tric polarizability from the lattice, the background Aeld 
method is employed, as described in Ref. ll^ l35l . The 
essence of this method is to measure baryon two-point cor¬ 
relation functions in the presence of a uniform electric Aeld 
£. Working in Euclidean space, the background held in¬ 
duces a quadratic Stark shift in the mass of the SU(4) 
ground-state baryon, 

Eb,4c = ± 2Cf\S\'^ ± O ( 1 ?'^) , (4) 

where Cf is the desired polarizabilitjj^ as deAned in 
Eq. Q. 

Due to the scalar nature of the SU(4) baryon ground 
state, this relation is equivalent to what one would expect 
for mesons. Eor comparison we also study the fermionic 
SU(3) baryon, whose energy shift contains an additional 
contribution from the non-zero magnetic moment /is ll35l . 

Eb,3c = ± ^2C'f — ^^3 ^ \£\ ± O . (5) 

Eor the SU(3) theory, we must therefore determine as 
well in order to extract Cf from the background held de¬ 
pendence. 


^ The electric polarizability of the neutron is more commonly defined in 
terms of the induced dipole moment d = giving a quadratic Stark 

shift of AE’n = • S = 2'KaE\S\^. In our notation a.^ = CfItt. 
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The background field method is implemented following 
Refs, ll^ [35]l where the uniform background field is in¬ 
cluded by multiplying the unitary gauge links by two phase 
terms, chosen so that the field is oriented in the z direction. 
Quantization of the uniform background field on a torus 
restricts the available field strengths to values 
, 27ren 

^ kiowliv^iv.’ ^ ^ 

where a is the lattice spacing, e is the electromagnetic cou¬ 
pling, Ng and Nt are the number of spatial and temporal 
lattice sites respectively, and giow is the lowest common 
denominator of the charges (for SU(3), qiov, = 1/3; for 
SU(4), glow = 1/2). 

For convenience we define a rescaled, dimensionless 
background field by = {ea^)~^£. We will analyze 
our lattice results using Eqs. Q-Q with all quantities re¬ 
placed with their rescaled, dimensionless counterparts, all 
of which will be denoted with a tilde: 


ms = ms/a, (1) 

Cp = Aircxci^Cp, ( 8 ) 

/is = i-iraflB- (9) 


Lattice details and fitting - The lattice calculations are 
done using the Chroma software package |[36l . We use the 
plaquette gauge action with unimproved Wilson fermions. 
The gauge configurations are quenched x Nt = 32^ x 
64 lattices (20000 heat-bath updates, 200 configurations 
separated evenly). For SU(4) we choose /3 = 11.028 
and for SU(3) /3 = 6.0175 following [371. Fermionic 
propagators are calculated for two different masses at each 
Npi value (tt = 0.1554,0.15625 for Np = 4 and k = 
0.1537,0.1547 for Np = 3), chosen such that the ratio 
of the pseudoscalar to vector meson masses mps/my = 
0.77 and 0.70 are matched llT4llT7l . 

Background field measurements are performed at six 
field values [n = 0,..., 5, see Eq. for both Np = 4 
and Np = 3, with correlation functions measured us¬ 
ing 40 evenly separated sources in (x, y) along the t = 
z = 0 plane. Each zero and non-zero field value has 8000 
measurements. All two-point correlation functions are fit 
over the range t G [4, 28] using fully correlated, multi¬ 
exponential fits including three excited states. 

Eor Np = 4, the two-point baryonic correlation func¬ 
tion in background field £ takes the form 


CBit,£) ~ Z{£)exp 


-tEB{£) 


( 10 ) 


at large t. Results for Eb{£) are then fit to Eq. (|^. We in¬ 
clude higher-order contributions from the background field 
following Il35l . 


Cf{£) = Cf + C'p\£\\ (11) 

The analysis for Np = 3 is complicated by the contribu¬ 
tion of the magnetic moment fip to the baryon self-energy. 
Eollowing ll35l . we make use of the boost projections 

'P± = ]^{1± i'js'ji ), ( 12 ) 
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FIG. 1. The ground state energy (in lattice units) vs. applied elec¬ 
tric field S for SU(4) baryons (top) and SU(3) baryons (middle), 
and ratio of projected SU(3) correlators vs. S (bottom). Their re¬ 
lations to the magnetic moment and polarizabilities are presented 
in Eqs. Q, (|^ and ( [T4| . Results shown are for the ensembles 
with mps/my = 0.70. 


and the boosted correlators 


cm 


{B(x,t)V±B{0,0))£ 
Z±{£)exp —tEpiZ) 


(13) 


The boost-projected amplitudes Z± contain equal and op¬ 
posite contributions from the magnetic moment, which we 
isolate by combining them in the ratio 


^ ZFi£)-Z_i£) 


£fiB{£) 

2m?p 


(14) 


Zf{£) + Z_{£) 

A simultaneous fit of Eb to Eq. ^ and the amplitude ra¬ 
tio in Eq. (14i allows us to determine both Cp and fip- 
To extract the polarizability Cp we use a fully correlated 
quadratic fit following Eqs. M- Once again we incor¬ 
porate quadratic terms to both Cp [as in Eq. (Ill] and fip, 

Ab(^) = As +(15) 


The polarizability results for both SU(4) and SU(3) are 
presented in Table [I] and results for the energies and the 
ratio Z^ vs. background field for the mps/rriy = 0.70 
ensembles are plotted in Eig. [T] 
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N]j 

mpsl'^v 

ihB 

OlC F 


(j-B 

Mb 

X^/dof 

4 

0.77 

0.98204(93) 

0.1420(56) 

-0.089(29) 

— 

— 

0.7/3 


0.70 

0.88805(113) 

0.1514(106) 

-0.142(68) 

— 

— 

4.8/3 

3 

0.77 

0.69812(51) 

0.2829(127) 

-0.177(45) 

-6.87(26) 

714(103) 

3.0/7 


0.70 

0.61904(59) 

0.2829(81) 

-0.165(24) 

-5.55(18) 

396(78) 

13.4/7 


TABLE I. Results for the polarizabilities and magnetic moments 
of the baryonic composites of a strongly-coupled SViNjj) theory, 
in lattice units. 


Constructing the dimensionless product aCprh^ (as 
needed for the DM cross section), we find that the SU(4) 
polarizabilities are larger than SU(3) hy about 50%. Thus, 
we find the SU(3) and SU(4) polarizabilities to be compa¬ 
rable when normalized to the baryon mass. Of course, the 
baryon mass itself scales proportional to N^', if we were 
to set the scale using a quantity such as the string tension 
which does not scale with Njj, then the SU(3) polarizabil¬ 
ity would be larger. 

The effect of the quenched approximation, in which dy¬ 
namical fermion loops are omitted from the lattice calcu¬ 
lation, is not straightforward to estimate. However, the ef¬ 
fects of such loops are expected to be suppressed with large 
Njj and heavy fermion mass; we note that even for QCD 
with its much lighter fermions, the effects of quenching are 
generally at most of order 10% ll38l . 

Our calculations are performed at a single lattice spac¬ 
ing and volume, both of which can lead to additional 
systematic effects. We expect all of these corrections to 
be small relative to the order of magnitude uncertainty 
taken for the nuclear matrix element Mp. As a cross¬ 
check, we note that the neutron polarizability from the 
PDG ||39l gives Cpm^ ~ 0.36 at the QCD physical ra¬ 
tio mps/my = 0.18, while our SU(3) lattice simulations 
give Cpm\ ~ 0.84 at mps/my = 0.70. These results 
are broadly consistent with the expected scaling of the po¬ 
larizability and baryon mass with mps- 

Direct detection cross sections — To relate the dimen¬ 
sionless lattice results to the dimensionful DM mass, mp, 
that we vary continuously in order to scan the parameter 
space of the theory, it is most convenient to give units to 
the lattice spacing a = ms/ms. Along with Eq. (j^, this 
leads to the physical value of the polarizability 

Cp = dvro: f Cp . (16) 

\mBj 


Putting everything together, the spin-independent cross 
section written as the conventional per nucleon rate for a 
nucleus with {Z, A) from Eq. (|^ becomes 


nucleon ( 


lUTra^plpiM^r 


A2 


m%R‘^ 


am^CpY , 
(17) 


where we use our lattice results in Table U to evalate the 
factor in square brackets. We emphasize that, unlike Higgs 
exchange, the cross section per nucleus scales as and 


not A^, and so the cross section per nucleon must be calcu¬ 
lated for each nucleus separately in order to compare with 
experiment. The strongest bound on the spin-independent 
direct detection scattering rate is from EUX |[T1. In Eig.[^ 
we show the scattering cross section per nucleon for xenon, 
and compare with the EUX bounds. We plot only the 
Np, = 4 case here, as the Np = 3 baryons are already 
excluded up to ~ 20 TeV in mass by the EUX bounds 
through their magnetic moments |[T^ . 

Discussion - Our lattice results have allowed us to 
calculate the spin-independent scattering cross section of 
SU(4) stealth DM through polarizability, which we com¬ 
pare against the EUX constraints in Eig. We find DM 
masses less than about 200 GeV are excluded, while the 
DM mass range 200-700 GeV could be probed by fu¬ 
ture experiments before reaching the neutrino background 
001. Currently, the strongest lower bound on the DM 
mass arises indirectly from the constraints on the lighter 
electrically-charged mesons that can be produced and de¬ 
cay promptly in collider experiments. Using our results 
ll23l , we estimate that DM masses below about 280 GeV 
are excluded given the EEP II bounds on the pseudoscalar 
mesons. 

It is remarkable that a composite DM particle with a 
weak-scale mass, composed of dark fermions charged un¬ 
der the weak and electromagnetic interaction, can never¬ 
theless be safe from both direct detection constraints and 
the EEP II constraint once mp ^ 300 GeV. This sug¬ 
gests there is a serious opportunity for future direct detec¬ 
tion experiments to probe the model. Given that the scat¬ 
tering cross section per nucleon scales as Z^jR'^) in 
Eq. ( [T7] l, the experiments with the heaviest nuclei are often 
more sensitive, i.e., xenon is 3.4 times more sensitive than 
argon if both experiments reach the same limit on the (con¬ 
ventional) spin-independent scattering per nucleon through 
Higgs exchange. 

With our lattice calculation of the dark matter polariz¬ 
ability in this model, the dominant remaining uncertainty 
stems from the treatment of the non-perturbative nuclear 
matrix element in Eq. Q, which is similar to the matrix el¬ 
ements required for double beta decay. A significant source 
of uncertainty is, for example, the presence of excited states 
in Xe-129 and Ge-73 that have energies of 30 and 15 keV, 
which will be probed by the loop in the cross section calcu¬ 
lation (typical momenta exchanges are roughly at the MeV 
scale). These resonances could appreciably change the re¬ 
sulting cross section, though the steep dependence on the 
dark matter mass suggests only a modest equivalent shift 
of the DM mass. 

The brightest opportunity for stealth dark matter discov¬ 
ery may fall within the domain of the Earge Hadron Col¬ 
lider (and future colliders). Meson phenomenology is very 
promising, since charged mesons can be produced through 
electroweak processes and decay completely into SM par¬ 
ticles. In contrast, production of the dark matter baryon is 
rare, since it is considerably heavier than the mesons and 
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FIG. 2. The DM spin-independent scattering cross section per nu¬ 
cleon evaluated for xenon is shown as the purple band obtained 
from the SU(4) polarizability, where the width of the band cor¬ 
responds to 1/3 < Mp < 3 from low to high. The blue curve 
and the light blue region above it is excluded by the LUX con¬ 
straints (il. The vertical, darker shaded region is excluded by 
the LEP II bound on charged mesons ||23]| . The orange region 
represents the limit at which direct detection experiments will 
be unable to discriminate DM events from coherent neutrino re¬ 
coil da. We emphasize that this plot is applicable for xenon, and 
would require calculating Eq. (13 to apply to other nuclei. 

would have form factor suppression. This implies the stan¬ 
dard missing energy signals that arise from DM production 
and escape from the detector are rare. 

Finally, there are many avenues for further investiga¬ 
tion of stealth dark matter, detailed in Il2^ . One vital is¬ 
sue is to better estimate the abundance. In the DM mass 
regime where stealth DM is detectable at direct detection 
experiments, the abundance of stealth dark matter can arise 
naturally from an asymmetric production mechanism Il23l 
that was considered long ago (THU and more recently re¬ 
viewed in BTI . If there is indeed an asymmetric abundance 
of bosonic dark matter, there are additional astrophysical 
consequences ll42ll44l that warrant further investigation to 
constrain or probe stealth DM. 
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